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oo ■ Abstract 



The construction of the multi-focal Veselago lens predicted earlier is proposed on the basis of a 
uniaxial photonic crystal consisting of cylindrical air holes in silicon that make a triangular lattice 
in a plane perpendicular to the axis of the crystal. The object and image are in air. The period 
of the crystal should be 0.44/zm to work at the wavelength 1.5/im. The lens does not provide 
superlensing but the half-width of the image is 0.5A. The lens is shown to have wave guiding 
properties depending on the substrate material. 

PACS numbers: 78.20.Ci,41.20.Jb, 42.25.-p 
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1 .INTRODUCTION 

A left handed medium (LHM), defined by Veselago Jl] as a medium with simultaneously 
negative fi and e, has recently attracted much attention mostly because of negative refraction 
at its interface with a regular medium (RM). This effect allows creation of a unique device 
called the Veselago lens. This lens is a slab of a LHM embedded inside a RM with the 
condition that both media have the same isotropic refractive index and the same impedance. 

It was shown recently^, [J that a two-dimensional (2D) dielectric uniaxial photonic crys- 
tal (PC) made of non-magnetic materials can behave as a LHM with negative e and /i if 
it has a negative group velocity in the vicinity of the r-point of the second Brillouin zone. 
Experimental demonstration of negative refraction in a metallic PC using the modes near 
the r-point has been presented in Ref. 0]. 

Proximity to the r-point is crucial for the PC-based LHM. Ingeneral, e and fi of a PC 



arc functions not only of u.\ but also of k. It has boon si i own (J that . in tin* vicinity of 
the r-point, e, ft, and uj are functions of k. Therefore, for every propagating mode, k, e, 
and fi can be represented as a function of uj. However, the parameters e(u) and fi{oj) are a 
property of a given mode in the medium rather than a property of the whole medium. It is 
shownP, 0| that parameters e and fi for evanescent waves (EW's) in the same medium and 
at the same frequency strongly depend on k. 

We think now that the 2D PC with a square lattice considered in Ref. 0,0,0] is a wrong 
choice from a practical point of view. Although both fi and refractive index n are zero[2j 
exactly at the r-point, one should still work near the r-point where the parameters e, fi, and 
uj are isotropic. In the square lattice the maximum value of n that can be achieved in the 
isotropic region is about 0.33. For the Veselago lens, the LHM and the RM should have the 
same value of n and it would be difficult to find a transparent RM with such a low n. The 
situation is much better for the 2D lattice with hexagonal symmetry in the plane (triangular 
lattice shown in Fig. 1(a)) because the width of the isotropic region near the r-point is much 
wider. Wang et al.y\ were the first to find out computationally that n = 1 can be reached 
in the isotropic region of the hexagonal (triangular) lattice. Their result can be explained 
from symmetry. For both square and triangular lattices, the expansions of frequency near 
the r-point have the same form up to the quadratic term, namely, uj 2 = uj" 2 — ak 2 . In this 
approximation, the frequency of the mode is isotropic. The question is, what is the order of 
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FIG. 1: (a) Triangular lattice of circular cylindrical holes in a dielectric matrix with e m = 12, 
(jL m = 1. The radius of the holes R = 0.35d, where d is the period, (b) The lattice with the basis 
and parameters as in Ref. 8. The radii of large and small air holes are R = 0.4<i and R = 0.13d, 
respectively. The dielectric matrix has e m = 12.96, \i m = 1. 

magnitude of the first anisotropic term in this expansion? Group theory provides a simple 
answer. In a square lattice with the symmetry axis C4, the largest anisotropic term is of 
the order of k 4 and has a form h^ky. In a triangular lattice with axis C§, this term is of the 
order of fc 6 and has a more complicated structure. Therefore the isotropic region is much 
wider in the latter case. Common sense also suggests that the anisotropy should be less in 
the lattice with six nearest neighbors than in the lattice with four neighbors. 

With this finding, the creation of the Veselago lens based upon photonic crystal at the 
frequency of long-distance communication (1.5/im) becomes a reality. However, even if we 
are able to match refractive index with air, effective |e| and \y\ of the PC would still be 
different from those in air. In this case, the impedances are not matched, reflection is not 
zero, and one gets a multi- focal lens considered in Ref. jjj. The advantage of the multi- focal 
lens is that its object domain is not limited to a distance L from the slab, as in the case of 
a matched Veselago lens. Instead, it is just a half-space in front of the infinite LHM slab. 
The disadvantage is that part of the energy flux gets lost due to reflection. 

In this paper we simulate the imaging of a point source by a uniaxial PC with the 
triangular lattice in a plane perpendicular to the axis of the PC. We demonstrate the focus 
in the case a > L, where a is the distance from the point source to the lens. In the case 
a < L, we find two foci behind the slab that are located in predicted places. 

Wang et al. also considered theoretically focusing in a similar slab matched with 
air. They claim observation of "superlensing" or even "unrestricted superlensing." We have 
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found "regular" lensing rather than "superlensing" and explain it by decay of the evanescent 
waves (EW's) in the PC. The same authors^ found out that by adding a basis to the lattice 
one can improve "superlensing" in the far- field regime. The best choice for improvement is 
shown in Fig. 1(b). This improvement strongly contradicts to the theory of the PC-based 
LHMjsiy]. This theory predicts that in the far-field region the intensity of the image for the 
infinite lens depends on the wavelength but, since the wavelength is larger than the lattice 
period d, it does not depend on the crystalline structure. We checked the statement of Wang 
et a/. [9] and have not found any improvement of the image for the PC with basis. 

Our paper is organized as follows. In section 2 we present all computational results 
concerning the image of the Veselago lens based upon the PC with a triangular structure 
with and without basis and compare these results with a previous theory of regular lensing. 
In section 3 we carefully study the behavior of the EWs in this PC. In sections 2-3 we 
assume that the PC is infinite in the direction of its axis. This is not realistic since to work 
at 1.5/im we need a PC with period about 0.4/rni. Even with current silicon technology 
one should not expect to get deep holes at such a scale. Therefore in section 4 we consider 
the propagation of waves through the PC slab with the depth of the holes less than the 
wavelength. We find that the slab exhibits a strong wave guiding effect at the working 
frequency and propose an explanation for it. 



2. IMAGES OF THE MULTI-FOCAL VESELAGO LENS 

In sections 2,3 we use for our simulation a uniaxial PC infinite in the direction of the axis 
of the crystal (z-axis) . Thus we have a 2D problem in the x — y plane. The cross-section of 
the PC by this plane is shown in Fig. 1. To begin, we discuss the lattice without basis (Fig. 
1(a)). All our simulations were built using the finite-element software package FEMLAB. 
Fig. 2 shows the spectrum of s-polarized propagating waves (H z = 0) for this PC. 

The dashed line shows our working frequency uj w d/2^c = 0.292, where d is the period 
of the crystal. At this frequency the refractive index n, as defined by equation n = ck/uu, 
is 1 and the wavelength A = c//0.292. Therefore one can get wavelength A = 1.5/xm under 
the condition that the period d = 0.44/xm and hole radius is 0.35c/. We have checked that 
the spectrum is isotropic at this frequency, which means that the equation u>(k x , k y ) = oj w 
is fulfilled on a curve very close to a circle in k x , k y space. Since the group velocity at this 
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FIG. 2: Six lowest bands of the photonic spectrum of the PC that has the geometry of Fig. 1(a). 

frequency is negative, it follows from Ref. Q] that both e and \i should also be isotropic 
and negative. However, they should not be equal to —1. The physical reasons for the 
appearance of a magnetization in the PC made of non-magnetic materials are discussed in 
Ref. |2|. Thus, at the working frequency, this PC is a LHM. 

We present a computer simulation of a lens made of the slab of uniaxial PC with a planar 
(x, y) structure shown in Fig. 1(a). The slab is surrounded by a homogeneous medium with 
e = fi = 1. It is infinite in the z-direction, has a width L = 3.83A in the x-direction and 
a length H = 17.52A in the ^/-direction. A 2D point source is located on the x-axis at a 
distance a in front of the slab. Thus, the problem is two-dimensional. 

Our results are shown in Figs. 3 and 4. In both figures part (a) shows the distribution of 
electric energy in the x — y plane behind the lens while part (b) shows the same distribution 
SIS db function of x/X at y — 0. 

The following theoretical results are necessary to interpret these data. It is shown in 
Ref. 0| that the positions of the foci behind the LHM slab are described by the equation 
L(2m — 1) — a with the origin at the right surface of the slab. In the case of a < L, 
m = 1,2, 3, ... and in the case of a > L, m = 2, 3, 4, .... The intensity of light in any focus is 
given by the equation 

I m = (1 - r 2 ) 2 r 4m - 4 / . (1) 
Here Iq is the intensity of the source and r is the reflection coefficient at a single RM-LHM 
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FIG. 3: Distribution of electric energy behind the PC slab for a = 0.5L in the x,y plane (a) and 
along the line y = (b). Two thin lines show the theoretical position of the two foci. Small arrows 
show computational positions of the foci. 







interface. If n = 1, the reflection coefficient is related to the electric permittivity as 

|e'| + e 1 ; 

and is independent of incident angle. It is important to mention that the intensity of the 
second focus at a < L is the same as the intensity of the first focus at a > L. Theory also 
predicts the appearance of foci in front of the slab that are created by reflected light, but 
we do not distinguish them because of the strong intensity of the source. 

The electric field in the ^-direction produced by the point source is E s = 
iE HQ 1 \pko) exp — ioot, where ko — w/c, n — 1, p is a polar coordinate with the origin 



at the source, and the Hankel function H 



(i) 



Jo + iN . The Fourier transform of the source 
field contains both the propagating and evanescent waves [5]. In the next section we show 
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FIG. 4: Same as in Fig. 4 for a = 1.5L. In this case we get only one focus. The second maximum 
is interpreted as a satellite of the focus. 

that the PC under study does not provide any amplification of EW's. Then we can consider 
only propagating waves in the Fourier transform to get the following expression |6] for the 
field near any focal point m: 



E m (x',y) 



£ rk H exp i (^ky + x' \Jk$ — k 2 — cut 
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dk. 



(3) 



The amplitudes E m < Eq take into account reflection and they can be related to the 
intensities I m given by Eq. (0) while x' is the x-coordinate with the origin at a focal point 
m. The wave vector ku — > k when the length of the PC slab H — > oo. For the finite lens 

(see Ref. |5|). The energy distribution given by \E m (x',y)\ 2 
provides the diffraction limit for the Veselago lens. 

To interpret the computational results shown in Figs. 3 and 4, one should take into 
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FIG. 5: Solid line shows the computational distribution near the first focus for the case of a > L 
and dotted line shows the computational distribution near the second focus for the case of a < L. 
There is no any numerical fitting in this plot. 

account that the field near the foci as given by Eq. (JHJ) has a main maximum at the focal 
point and infinite number of smaller maxima. It is also important to notice that due to 
reflections the amplitude E m strongly decreases with increasing m. Fig. 3 (a < L) shows 
the first maximum (m = 1) exactly at the theoretical position and a long plateau with two 
small maxima. We think that the first of these two is a satellite of the strong maximum 
with m = 1 while the second one corresponds to m — 2. It has a small shift to the left from 
its theoretical position. 

In the case a > L, the matched Veselago lens does not produce any foci. The foci appear 
due to reflection and the first focus should be at a distance 3L — a from the lens. The first 
maximum in Fig. 4 is very close to the theoretical position and we think it is the focus. We 
interpret the second maximum in Fig. 4 as a satellite of this focus. 

The following arguments confirm such an interpretation for a < L and a > L: 

a) The intensity of the focus at a > L is very close to the intensity of the second focus at 
a < L. Eq. (JTJ) predicts that the intensities should be equal because both additional foci 
are formed by equal amount of reflections. 

b) Fig. 5 shows the lateral (^-direction) distribution of energy for both additional foci. One 
can see that the main maxima have almost the same shape. Note that no fitting parameters 
are used in Fig. 5. 
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FIG. 6: Distribution of electric energy along the lateral direction near the focus for the case of 
a > L. Solid line shows the computational distribution near the first focus and dashed line show 
the analytical result for intensity obtained from Eq. (J3J). The results are normalized to have the 
same maximum value at the focal point. 
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FIG. 7: Distribution of electric energy along the lateral direction near the first focus for the case of 
a < L. Solid line shows the computational distribution for PC as shown in Fig. ^a), dash-dotted 
line shows the computational distribution for PC with "basis" as shown in Fig. ^b) and the dashed 
line shows the analytical result for intensity obtained from Eq. ©. The results are normalized to 
have the same maximum value at the focal point for PC without "basis" . 
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c) In Fig. 6, the computational distribution of energy near the additional focus at a > L in 
the lateral direction is compared with the analytical expression given by Eq. (jSJ) • The value 
of E m at m = 2 is chosen in such a way that the maximum values of both distributions 
are the same. The main maximum of the distribution is described satisfactorily by this 
equation. 

d) We believe that the second maxima in both Fig. 3(b) and Fig. 4(b) are satellites of 
the previous foci rather than additional foci because they are located at the same distances 
from the previous foci. The minima on both curves are also located at the same distances 
from the previous foci. This similarity indicates that the first and the second maxima are 
described by the same universal curve, as it should be if the second maximum is a satellite 
of the first one. 

Now we describe the intensity distribution in the vicinity of the main focus of the Veselago 
lens that exists at a < L only. The lateral distribution of energy is shown in Fig. 7. One 
can see that the diffraction limit given by Eq. (|3~|) is a little narrower than the peak obtained 
by simulation. Thus, we have regular lensing rather than superlensing. The main peak 
obtained on the lattice with a basis is indistinguishable from the result without a basis. 

The simulation gives much sharper peaks in the lateral direction than in the perpendicular 
direction. This is a general property of Eq. Q that has been observed in the case of the 
completely matched lens as well|5|, |6|. Note also that in the lateral direction the simulation 
does not reproduce the satellites. 

We think that the difference between computational and analytical distributions of the 
satellite maxima in the lateral direction that we see in Figs. 5, 6, and 7 is due to some 
angular dependence of the reflection coefficient r that is absent in the theory (see Eq. (j2J)). 
This dependence may appear because our working frequency is not close enough to the T- 
point and the PC cannot be perfectly described as a continuous medium. It would certainly 
affect the fine structure of the image more strongly. Note that the simulation describes the 
satellites much better in the case when reflection is smally, Igl. 



3. EVANESCENT WAVES IN TRIANGULAR PC 

The initial idea of superlensing was proposed by Pendry[ll] It is based upon amplifica- 
tion of evanescent waves by the LHM slab that should help to restore the image. Following 
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VeselagoQ, Pendry considered a hypothetical LHM (HLHM), which is a homogeneous ma- 
terial with simultaneously negative e and u. This amplification is shown to be related to 
surface polaritons[3| studied by Ruppin 13] in the HLHM. The polaritons are closely con- 
nected with negative e and \x in the bulk of the HLHM. 

It has been shown recently that a two-dimensional PC with a square lattice generally 
does not amplify EW's. However, some amplification in the near-field regime may or may 
not appear depending on the properties of the surface. For example, if the surface is cut 
between the holes (BH), amplification does not appear at all, while if it is cut across the 
holes (AH), amplification appears and provides superlensing in the near- field regime. The 
authors argue that the amplification is due to surface waves that has nothing in common 
with the left-handed properties of the PC with respect to propagating waves. 

We show here that a 2D PC with the triangular lattice does not amplify EW's for both 
BH and AH cuts of the surface. The geometry of the simulation is shown in Fig. 8 and 
the working frequency is the same as in the previous section. The PC slab is surrounded 
by a homogeneous medium with e — fi — 1. An EW in a homogeneous medium has a form 
e z ~ exp(ik y y — kx) with the dispersion law k 2 — k 2 = u 2 /c 2 . We create such a wave in our 
simulation using boundary conditions for an incident wave at x = in a form e z = exp (ik y y) 
and de z /dx = —ne z . Note that due to reflection the total field shown in Fig. 8 may not 
coincide at x = with the field of the incident wave. The coordinate of the left surface of 
the PC slab is x = L/2. 

The result shows that the EW decays inside the PC slab without any signature of ampli- 
fication. The surface of the PC slab is BH. We have checked that the PC slab with the AH 
surface gives the same results. We have also performed similar simulations for the PC with 
the basis shown in Fig. 1(b) with the same result. 



4. WAVEGUIDE EFFECT IN THE LENS WITH A SMALL SIZE ALONG THE PC 
AXIS 

We have considered so far only PC slabs that are infinite in the axial z-direction. In 
this section, we investigate the effects of termination of the slab in the z-direction. Namely, 
we are interested in how power in the x-direction varies as a function of the distance from 
the slab interface. We still keep our crystal infinite in the y-direction by using periodic 
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FIG. 8: Electric field of the evanescent wave. The cross-section y = is shown. The solid line is 
for k = O.lko and k y = v/LOl^o, the dotted line is for k = 0.2fco and k y = VLMko, and the dashed 
line is for k = ko and k y = \/2&o- For other k in the interval 0.2ko and ko, the plots are between 
the dotted line and dashed line. The surface cut is BH. The plot for AH cut is indistinguishable 
from this plot. 

boundary conditions. We are not interested now in the losses due to reflection of waves from 
the interface. Therefore we define the power P as the surface integral 

P{x) = [ Sda, (4) 
J Q 

where S is Poynting's vector in the x-direction and the surface Q is perpendicular to this 
direction and located at a distance x from the left surface of the PC slab (see Fig. 9(b)). 
Calculating this integral as a function of x inside the PC per unit length in ^/-direction one 
can study the leaking of energy due to diffraction in the z-direction. Note that we ignore 
small absorption of energy in pure silicon so that diffraction is the only mechanism that 
provides the x-dependence of P. The diffraction depends on what kind of material is above 
and below the PC slab in the ^-direction. One side of the slab is always surrounded by air 
while the other side adjoins different media, namely air, silica or silicon. 

A 2D, top- view representation of our model is shown in the x — y plane in Fig. 9(a). 
Periodic boundary conditions in the y-direction were used to simulate an infinite slab in 
this direction. As depicted in Fig. 9(b), an s-wave is incident normally through the slit 
in the boundary with low- reflecting boundary conditions (LRBC). The incident wave then 
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FIG. 9: Representation of the PC slab geometry, (a) Top view of the cross-section of the PC in the 
x — y plane. The system is infinite in y-direction due to periodic boundary conditions, (b) Side- 
view. Normal light is incident through a slit on the left while the remaining external boundaries 
are low-reflecting. Normal power flux is calculated across the surface Q indicated by dashed line. 

goes through the air buffer to the PC slab surface. The transmitted and reflected field is 
almost perfectly absorbed by the surface with LRBC surrounding the system. We measure 
the power transmission in the x-direction as a function of distance in the slab. The size of 
our model is limited by the computational memory requirements for large 3D geometries; 
we can accurately model a slab about seven lattice cells thick. We use AH cut of the surface 
with length in the x-direction 3.54A and with height in the z-direction 0.96A. 

The results of our simulation are shown in Fig. 10. One can see that the slab surrounded 
by air on all sides is a very good waveguide. The loss of electromagnetic power leaking from 
the slab in the z-direction of about 1.27% per wavelength or 0.38% per lattice cell (-\/3d). 
If, however, the slab were replaced by air, power flux transmission drops significantly faster 
due to regular diffraction from the slit. Thus the PC slab suppresses diffraction. We think 
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FIG. 10: P(x)/P(0) in percent in the PC surrounded by air (solid line circles) and on substrates 
of silicon (solid line triangles) and silicon dioxide (solid line squares). The dash-dotted line shows 
the same power ratio when the PC slab is substituted by air and the "substrate" is also air. The 
trend line (dashed line) indicates least squares linear fit of the PC in air data. 

that the suppression is due to total internal reflection while all remaining losses are due to 
dipole radiation from the PC surface. A very low refractive index n — 1 for the propagating 
waves is only because of the proximity of the wave vector with k z = to the T-point. A 
two-dimensional band diagram shown in Fig. 2 does not have any meaning for a diffracted 
wave with k z ^ 0. The diffracted waves, therefore, should have a much larger value of n and 
should experience total internal reflection at the interface with air. To check this hypothesis 
we simulated the PC slab on a substrate made from silica or silicon. The silica substrate 
(n = 1.45) does not prevent wave guiding, but the silicon substrate (n = 3.6) makes it 
worse. The losses with the Si-substrate become about 1.82% per lattice cell. Power loss, 
however, is still much smaller than one would expect from diffraction; perhaps there is an 
additional mechanism for wave guiding. Thus, both substrates are good candidates for the 
construction of a Veselago lens. 



5. CONCLUSION 



To summarize, we propose to make a Veselago lens using a uniaxial PC with a triangular 
lattice of cylindrical holes as a LHM. The working frequency is chosen in such a way that 
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the refractive index for propagating modes n = 1. Using air as a regular material one can 
match n, but cannot match the impedance. The lens, therefore, is multi-focal. The EW's 
decay inside the PC and super lensing is absent. However, the half- width of the image is 
about 0.5A, in agreement with the diffraction limit for the flat lens. The additional foci are 
more smeared. Unexpectedly, the lens has good wave guiding properties that are mainly 
due to total internal reflection, the quality of which depends on the substrate. 
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